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Abstract 

We study a new simple quadrature rule based on integrating a quadratic 
spline quasi-interpolant on a bounded interval. We give nodes and weights for 
uniform and non-uniform partitions. We also give error estimates for smooth 
functions and we compare this formula with Simpson’s rule. 


1 Introduction 

We study a new simple quadrature formula (QF) based on integrating a quadratic 
spline quasi-interpolant (qi) on a bounded interval. It can be considered as the 
companion of Simpson’s rule (as midpoint rule is the companion of trapezoidal rule) 
in the sense that the signs of errors for both QF are opposite in almost all examples. 
Here is an outline of the paper: in section 2, we recall the dehnition and main 
properties of the quadratic spline quasi-interpolants which are used. In section 3, we 
give the weights of the QF. In sections 4 and 5, we give estimates of the quadrature 
error for smooth functions in the cases of non-uniform and uniform partitions. In 
the latter case, we develop a more deeper study by using the associated Peano kernel 
and we compare our results to those of Simpson’s rule. Finally, we illustrate our 
results with two numerical examples. 

Regarding some recent papers close to ours, e.g. De Swardt & De Villiers jl], and 
Lampret nn, we use quadratic splines as the first one, however we do not consider 
interpolants, but quasi-interpolants (other QF based on spline quasi-interpolants are 
given in nHi) and our QF seems to have slightly better properties, at least than the 
Lacroix rule. The paper by Lampret is more oriented to the calculation of sums via 
the Euler-Maclaurin formula: it would be interesting to compare the performances 
our formula in this domain. 


2 Univariate quadratic splines and discrete 
quasi-interpolants 


Let X = Xn = {xq,Xi, ... ,Xn} be a partition of a bounded interval I = [a,b] , 
with xq = a and Xn = b. For 1 < i < n, let = Xi — Xi^i be the length of the 
subinterval Jj = Xi] and let F = F„ = {0,1,..., n + 1}. We denote by 52(X) 
the n + 2-diniensional space of quadratic splines on the partition X. A basis 
of this space is formed by the family of quadratic B-splines B = {Sj, z G F}, with 
triple knots a = a :_2 = x_i = Xq and b = Xn = Xn+i = Xn+ 2 - With these notations, 
the support of Bi is [a:i_ 2 , for all z G F. Define the set of data sites (or Greville 
points): 

©n = {Oi = +^i) D ^ r}. 

Note that 6*o = Xq and 9n+i = Xn- In jUjlHI, we have proved the existence of a 
unique discrete quasi-interpolant (abbrev. dQI) of type 

n+l 

i=0 

whose discrete coefficient functionals are respectively poif) = f{xo), /Un+i(f) = 
f(xn) and, for 1 < z < zz : 

hi(/) = + bif{9i) + Cif{9i+i), 


which is exact on the space 112 of quadratic polynomials. Using the B-spline expan¬ 
sion of monomials and setting efc(x) = x^ for fc > 0, we get (see e.g. |T].[T^.pnjl 

Co = ^ Bi, ei = ^ 9iBi, 62 = ^ Xi_iXiBi, 

ier ier zer 

and writing the equations Qck = Ck for 0 < fc < 2, we obtain: 


n?. 


—, = 1 Ci = 




with 


(Ti + 

hi 

Oi = 




-, a'i = \-ai = 


Oi + 
hi—I 


hi-i -\- hi * hi-i -\- hi 

We can now define the fundamental functions of S 2 {Xn) associated with the quasi- 
interpolant Q: 

Bq = Bq + UiBi, Bn+l = CnBn + Bn+1, 

and, for 1 < z < zz: 

Bi = Ci-iBi-i + biBi -\- aj+iSj+i. 

They allow to express Qf in the following shorter form: 


zer 






and to show that the inhnity norm of Q is eqnal to the Chebyshev norm of its 
Lebesgue fnnction: 

1 ^* 1 - 

ier 

In [20], Marsden proved the existence of a nniqne Lagrange interpolant v G 52(X) 
satisfying v{6i) = f{0i) for 0 < z < n + 1. He also proved the following interesting 
resnlt: 

Theorem 1. The inhnity norm of the Lagrange operator is nniformly bonnded by 
2, for any partition X of I. 

There exists a similar resnlt for the dQI Q dehned above : 

Theorem 2. For any partition X of I, the inhnity norm of Q is nniformly bonnded 
by 3. If the partition is nniform, one has ||Q||oo = fof ~ 1-4734. 

proof : For the nniform case, see e.g. in|. The hrst part of the proof is easy 
becanse we have 

\Bi\ < (|ci_i| + bi + \ai+i\)Bi, 

therefore, by snmming on all indices, as B-splines snm to one : 

^ \Bi\ < max(|ai| + 6* + \ci\) = max(l + 2(|aj| + |ci|)) = max(l + 2cri(j') < 3, 
since 6* = 1 — a* — q = 1 + \ai\ + \ci\, and cTj, a' < 1 □ 

Remark. The resnlts of this section are also valid when X contains some knots of 
mnltiplicity 2 or 3. Assnme hrst that ^ = Xp = Xp+i is a double knot, then Qf is 
only continnons at that point. Moreover as hp+i = 0, we have supp{Bp) = [a:p_ 2 ,^], 
supp{Bp+i) = [^-hp,^ + hp+ 2 ], supp{Bp+ 2 ) = [C,a;p+ 3 ]- Similarly, as cTp+i = 0, we 
have Op+i = Cp+i = 0 and = 1, hence: 

p n+l 

Qf = + fiOBp+i + Y. 

2=0 i=p-\-2 

Now, if p = Xg^i = Xq = Xq+i is a triple knot, then Q 2 / has a discontinnity at 
this point. Assnme that / is itself discontinnons and admits left and right limits 
f{p~) and Then as hq = hg+i = 0, we have supp{Bq) = [p — hq_i,p], with 

Bq{p-) = 1 and Bq{p+) = 0, while supp{Bg+i) = [p,p + hg+ 2 ], with Bg+i{p-) = 0 
and Bq^i{p~^) = 1. As aq = cr^+i = 0, we get Ug = Og+i = Cq = Cg+i = 0 and 
bq = bq^i = 1, hence: 

g—l n+l 

Qf = + fiV~)Bg + f{p^)Bq^l + P^{f)B,. 

i=0s i=q-\-2 


Finally, from theorem 2 and standard arguments in approximation theory (see [H]) 
we deduce : 


Theorem 3. There exists a constant 0 < C < 1 such that for all / G and 

for all partitions of J, with h = maxhj, 

\\f-Qf\\oo<Ch^\\D^f\\^. 

3 Quadrature formula associated with Q 

In this section, we compute the weights of the quadrature formula (abbr. QF) 
associated with the dQI Q, in the general case (non-uniform partition of the interval) 
and in the interesting particular case of a uniform partition. 

3.1 Case of a non uniform partition 

The QF associated with Q is of course obtained by integrating Qf\ 

IqU,I)= fQf = f{xo) + + 

J a J a J a J a 

We know that Wi = Bi = |(hj_i -|- 4hj -|- hj+i) for 2 < i < n — 1 and 

Wo = Wi = + h2), Wn = + 3hn), Wn+l = 

do D 6 

Finally, we get the following QF: 

ier 

which can also be expressed in the classical form: 

ier 


with Wo = Wo + aitci, Wn+i = -1- Wn+i-, and for 1 < i < n : 


Wi= Bi = d-iWi-i + biWi + tti+iWi+i. 

J a 

The explicit expression of these weights is rather complicated. We only prove the 
following result : 

Theorem 4. 

1) For all partitions, the weights Wi satisfy Xlier — 3(& — a)- 


2) More specifically, if there exists a r > 0 such that - < 

r 

n — 1 , then we have the more precise upper bound : 


hi+l 

hi 


< r for all 1 < f < 


\wi\ <{b- a)(l + 2 (^-)^) 
r + 1 


proof : \wi\ = ^.gp \jBi\<J ^.gp \Bi\ < max(|ai| + bi + |q|) / ^.gp Bi < 

{b — a) max(|aj| + 6 * + |q|) < 3{b — a), the last inequality being a consequence of 
theorem 2. In addition, when ^ < r for all i, then we immediately deduce 

that ^ < cri,cr- < ^ and \ai\ + bi + \ci\ = l + 2 (|ai| + |ci|) = l + 2 aia- < 1 + 2 (^) 2 , 
whence the result □ 


Remark 1. As Q is exact on 112 , we have Egi^f, I) = X(/, I) — I) = 0 for all 
/ G 112 , therefore we can deduce that Eg^f, I) = 0{h?), with h = max/ij for smooth 
functions / (see section 5 below, theorem 5). However, when the partition X satisfy 
the symmetry property Xi = Xn-i, 0 < i < n (e.g. Chebyshev knots or, more gen¬ 
erally, zeros of classical orthogonal polynomials in J = [— 1 , 1 ], see e.g.p][7jjH][Tn]h 
then the weights are also symmetric and we also obtain Eg{e^, I) = 0. In that case, 
the error is a 0 {h‘^) for / G where h = max hi. 


Remark 2. In general, the weights Wi are positive. However, they can be negative 
when the partition is strongly non-uniform. Take for example the following sequence 
of knots : 

X = X 7 = {-1, -0.9, -0.3, -0.2, 0.5, 0.6, 0.95,1} 
whose associated sequence of steplengths is : 


hx = {0.1,0.6,0.1,0.7,0.1,0.35,0.05} 


Here we have r = 7, 
bound J 2 ier < f| 


1 h- 

i.e. - < —— < 7 for all i, and theorem 4 gives the upper 

7 hi 

~ 5.06. The sequence of weights is the following : 


w = (.0146, .0122, .7463,. - .0622, .8780, -.0257, .4287, .0007, .0074} 


We observe that the weights W 3 and are negative : however, their absolute values 
are small. Moreover, we have X^jgp = 2.17 < 5.06. 

3.2 Case of a uniform partition 

Assume that n > 5 (in order to avoid boundary effects) and that the partition 
is uniform: hi = h = — for 1 < f < n. In that case, the weights can be computed 
once for all and we obtain, by setting fj = f{9j),j G T (see e.g. [TT)] [ T7] [ TH] h 

1 7 73 

h) = h[ -(/o + fn+l) + g(/l + fn) + ^(/s + fn-l) + ^ fi]- 





As for some Newton-Cotes QF, this QF has a better degree of precision: 

Theorem 4. The QF associated with the dQI Q on a uniform partition of I is ex¬ 
act on the space Ila of cubic polynomials. Therefore E{f,I) = O(h^) for / smooth 
enough. 

proof : This is due to the symmetry of weights and nodes with respect to the 
midpoint of /. □ 


4 Error estimate on a non-uniform partition 

As we have already seen in section 4, Eq^J,!) = —XQ{f,I) = 0{h?) when 

the partition X of / = [a, b] is non-uniform, with h = maxi<j<„ h^. In this section, 
we give more specihc results when / G From theorem 3, we immediately 

deduce: 

Theorem 5. There exists a constant 0 < C 3 < 6 — a such that for all / G 
and for all partitions X of /, with h = maxi<j<„ hf. 

|£«(/,-r)l SCsA’IIbVIIoo. 

proof : we know (theorem 3) that there exists a constant 0 < C* < 1 such that 
11 / - <3/lloo < Cli^WD^fW^. As EQ(f,I) = I(f - Qf), we can write: 

\EQ(f,I)\ < [ \f-Qf\ < (b-a)Ch^\\D^f\U 

J a 

Therefore we get the result with = {b — a)C. □ 

5 Error estimate on a uniform partition 

Assume that the partition X„ is uniform = h = — for 1 < i < n) and that 
/ G As is exact on IIs, the Peano kernel theorem (se e.g. j2], 

chapter III) gives: 


where 


EQ(f,I)=^-j K{t)D'f{t)dt, 


K(t) 



Q{. - t)\]dx, 


Q(- - 

ier 


and 


As fl^{x — t)\dx = JI’{x — t)^dx = \{b — t)^, we obtain: 

K{t) = ^(6 - tf -^Wi{ei - t)l. 
ier 

We see immediately that K{a) = K{b) = 0. First, it is clear that K{b) =0 since 
(Oi — 6 )+ = 0 for all i G F. Second, let p{x) = (x — a)^, then T(p) = j(b — a)^ = 
^q(p) = sincep G IIs, therefore iF(a) = = 0. 


5.1 Sign structure of the Peano kernel 

For the sake of simplicity, we now assume that / = [0,1], then h = ^ and 9i = {i—\) h 
for I < i < n, with 0o = 0 and 9n+i = 1. Therefore we have: 


K{t) 


ii-ty 








3 

+ 


+ 


7 

8 



n—2 


+ Q(i-'’ki+E 


i=3 



It is not difficult to prove that K{1 — t) = K(t). Actually, setting 


K{t) 


n+l 


-{l-tf-'^Wi{9i-t) 


i=0 


3 

+ 


and using the symmetry of nodes and weights, we get 


, n+l 

^(1 -t) = - 9i)%. 

i=0 

Now, we observe that {9i —1)\ — {t — 9i)\ = {9i — t)^ and that the cubic polynomial 
p{s) = [s — ty is exactly integrated by the QF, hence: 


n+l 


n+l 


p{s)ds = '^Wip{9i) -[(1 - - ty - '^Wi{9i - tf = 0 . 


i=0 


i=0 


The above properties imply that 


n+l 


K{t) - K{1 -t) = -[(1 - tf - tf - ^Wi[{9i - tf_^_ - (t - 9iff = 0. 


i=0 


In the interval h = [0, 9i] = [0, |], we have: 






K[(t) = 



Therefore, it is clear that Ki has a minimum that Ki{^) = 0 and 

K,{e^) = i^id) = Moreover K[{\) = §. 


In the interval I 2 = [9i, 62 ] = [|, ^], we have: 


A'(i) = A-,(i) = A'.(«)-y h-tf , A-'(i) = A;(«) 


21 h 


h 

t - 

2 


As K[{ 92 ) = K[{^) = we see that -^ 2 (^ 2 ) = 0, therefore we can factorize: 

f 3/i\ / 2 35h 


7h^ 
16 


and we deduce that = 0 for t = ^ and t = t = h ^ l.OSQh. Therefore 

K 2 has a maximum at the latter point and 

K2{t) ^ 0.03h" < 

moreover, we have -^ 2 (^ 1 ) = Ki{9i) = -^ 2 (^ 2 ) = and we can factorize: 


Mt) 





In the interval = [ 6 * 2 , 6*3] = [^, ^], we have: 


Kit) = Ksit) 


K2it) 


73h 

~n 



Kit) 



As ^' 2 (^ 2 ) = 0, we see that ^' 3 (^ 2 ) = 0, therefore we can factorize: 


Kit) 



9h 


t + 5A 



it — 2 h) (t 





and we deduce that has one maximum and two minima: 


K3i2h) 




64 ■ 


Finally, we also obtain the factorization: 


Ksit) 


64 


1 

4 



Now, in the intervals 9i\ = [(i — |)/i, (i — |)/i], for 4 < i < n — 2, if we 

assume that 



as we have, by definition 


we immediately obtain 


Therefore, we deduce that, in the interval J^, K = Ki has one maximum and two 
minima 

A'i((i - l)ft) = A'j - |)ft) = A. ((* - i)ft) = 

In the three last subintervals In-i = [1 — ^, 1 — In = [1 ~ ~ |] and 

Jn+i = [1 — |, 1], the behaviour of K is symmetrical of that one in the three first 
subintervals. To sum up, we obtain 

Theorem 6. The Peano kernel K is negative in the two small subintervals Ji = 
[ 0 , and J 3 = [1 — 1 ] and positive in the subinterval J 2 = [^, 1 — ^]. 

5.2 Error estimate for the quadrature formula 

Let / G C'^(/) be a given function, then from 


Eq(/,/) = - Kit)D^fit)dt 


and theorem 6, we deduce: 


Theorem 7. For any function / G there exists two points c,c & I such that 

EqU,!) = - h^D^f(c) - h^D^fic). 

J 57gQ M ; ^g2 M ; 

proof : From the mean value theorem, we know that there exists 
Cl G Ji, C2 G J 2 and C3 G J3 such that 


4h 

Eo(f,0 = ^ D"f(ci) ^ ■ K(t)dt + A>‘/(c2) 


K{t)dt + Defies) / K{t)dt 


We compute successively 


K{t)dt = 


I . f Ah\ , 

4 /n \ 9 295245’ 



I K{t)dt = / \t——]dt = 


2 


K{t)dt = 


9 J 37791360’ 
59h^ 


and for 2 < i < n — 2\ 


2880’ 

23h^ 


/ K{t)dt = 
'(i-i)h 960 


Then, nsing nh = 1 and setting 71 = 2 ^ 1 ^, 72 = ^,73 = we obtain 

snccessively: 




, ,23h^ 23h^ 23h^ 

' K(t)dt = in — 3) -=-. 

ih ^ ^ ^ 960 960 320 


K{t)dt = 72^^ — [2 X + 2 X 7777^ — 7777 ) = 72^.^ — 73^^) 


' J 2 


37791360 


2880 320 


/ K{t)dt = / K(t)dt = —7ih®. 
hi Jja 

Now, setting 7 ^ = 7^/6 for 1 < j < 3, we have: 

So(/. /) = 1 f K(t)D*f(t)dt = iJi*D*f(c2)-h‘ [7i£>7(ci) + 73C7(c2) + liD'fic^)] 

Finally, setting 74 = 271 + 73, the mean-valne theorem implies that there exists a 
point C 4 G / snch that 

EqU, I) = - 74h'ZlV(c4), 

with 72 = ^ and 74 = ^. □ 

Remark. For n large enongh, one can write 

23 

EQ{f,I) = ^h^Dh{c) + 0{h^), 
therefore, the main part of the error is contained in the hrst term. 


5.3 Comparison with composite Simpson’s rule and extrap¬ 
olation 

Let ns compare the above error with the qnadratnre error Es{f, /) = X(/, I) — 
Xs{f,I) of composite Simpson’s rnle. Taking n = 2m even, the latter is based on 
the n+1 points {xj, 0 < j < n}, while onr QF is based on the set {6*,, 0 < i < n+1} 
which has n + 2 points. It is well known (see e.g. 12)00^2]) that for any fnnction 
/ G C^(/), there exists a point d G / snch that: 


Es{f,I) = 


180 


h^Dhid). 















Therefore, when D‘^f is of one sign over the interval I, the quadrature errors are of 
opposite signs, i.e. the two formulas give lower and upper estimates of the value of 
the integral. 

Remark: as shown by numerical examples below, the following linear combination 
obtained by extrapolation 

iQsif, I) = 4^32Jq(/, I) + 23Xsif, /)) 

gives a still better approximation, with an order 0{h^) (notice the nice symmetry 
of weights). 


6 Numerical results 

For the two following functions /i ,/2 and /a (other examples are given in jIHI), we 
give the quadrature errors Eq^J), Es{f) and EQs{f) in terms of the number n of 
subintervals: 


Example 1: 

X{fi, I) = /q IQx^/'^sin (a:2) dx = 3.2523064663781227544. 


Example 2: 


X(/2,/) 



n 

EqUi) 

EsUi) 

EqsUi) 

64 

-.86(-7) 

1.23(-7) 

1.13(-9) 

128 

-.54(-8) 

.76(-8) 

.16(-10) 

256 

-.34(-9) 

.47(-9) 

-.40(-12) 

512 

-.21(-10) 

.29(-10) 

-.52(-13) 

1024 

-.13(-11) 

.18(-11) 

-.33(-14) 

1 

1 

0.8 


0.3)2 

+ 0.01 ' (x- 

-0.7)2 + 0. 

04/ 

n 

EqU2) 

Es{f2) 

EQs{f2) 

64 

-.19(-5) 

.23(-5) 

-.14(-6) 

128 

-.ll(-6) 

.14(-6) 

-.37(-8) 

256 

-.67(-8) 

.90(-8) 

-.ll(-9) 

512 

-.41(-9) 

.56(-9) 

-.35(-ll) 

1024 

-.25(-10) 

.35(-10) 

-.11(-12) 




Example 3 : 


X(/3,/) = 


/_! 1 + 16x‘^ 


dx = 0.6629088318340162325296195. 


n 

EqU2) 

EsU2) 

EqsU^) 

256 - 

-.33(-10) 

.46(-10) 

-.44(-12) 

512 - 

-.21(-11) 

.28(-ll) 

-.13(-13) 

1024 - 

-.13(-12) 

.18(-12) 

-.42(-15) 

2048 - 

-.80(-14) 

.11(-13) 

-.13(-16) 

4096 - 

-.50(-15) 

.69(-15) 

-.41(-18) 

of errors 

are opposite when h 

is small enough. 
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